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Abstract 

We have used the locally self-consistent Green's function (LSGF) method 
in supercell calculations to establish the distribution of the net charges as- 
signed to the atomic spheres of the alloy components in metallic alloys with 
different compositions and degrees of order. This allows us to determine the 
Madelung potential energy of a random alloy in the single-site mean field 
approximation which makes the conventional single-site density-functional- 
theory coherent potential approximation (SS-DFT-CPA) method practi- 
cally identical to the supercell LSGF method with a single-site local inter- 
action zone that yields an exact solution of the DFT problem. We demon- 
strate that the basic mechanism which governs the charge distribution is the 
screening of the net charges of the alloy components that makes the direct 
Coulomb interactions short-ranged. In the atomic sphere approximation, 
this screening appears to be almost independent of the alloy composition, 
lattice spacing, and crystal structure. A formalism which allows a consis- 
tent treatment of the screened Coulomb interactions within the single-site 
mean-filed approximation is outlined. We also derive the contribution of 
the screened Coulomb interactions to the formalism and the general- 
ized perturbation method. 
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I. INTRODUCTION 



The coherent potential approximation (CPA) [|I]-§J as implemented on the basis of 
multiple scattering theory [|],[5] and combined with density functional theory (DFT) ||-§| 
constitutes the basis for ab inito calculations of the electronic structure and physical 
properties of random metallic alloys. This combination of the CPA with DFT, or, in 
most cases, with the local density approximation (LDA), seems to be quite transparent 
MM leading to expressions for the one-electron potential and total energy which are 
very similar to those for ordered systems. However, there is, by now, a well-recognized 
problem [pD|-[TB[] with this description caused by the fact that the atomic or "muffin-tin" 
spheres, which artificially divide the crystal into regions associated with particular alloy 
components, may poses non-zero net charges. 

The problem stems from the fact that the conventional single-site (SS) DFT-CPA 
method is based on the effective medium model of a random alloy which considers only 
conditionally averaged quantities and leads to the use of the single-site approximation 
not only in the electronic structure part of the problem during the solution of the CPA 
equations, but also in the DFT self-consistent loop in the calculations of the electrostatic 
contributions to the one-electron potential and energy. The single-site approximation pro- 
vides no information as to the charge distribution beyond the atomic sphere of each alloy 
component and, since the surrounding effective medium is electroneutral, Poisson's equa- 
tion cannot be solved properly if the atomic spheres have non-zero net charges. Hence, to 
find the correct solution to Poisson's equation one must somehow describe the effect of the 
missing charge. Since the electron density inside each atomic sphere is well-defined, any 
such description may be associated with a modification of the effective medium specifically 
for each alloy component. This may be regarded as an inconsistency since, in that case, 
the CPA and the electrostatic part of the DFT are based on different effective media. 

One obvious solution to the problem is to use electroneutral spheres (see, for instance, 
Ref. [|13|]). However, in the methods based on the atomic sphere approximation (ASA) 
this frequently leads to large sphere overlaps and a quite poor description of the electronic 
structure, especially in the case of inhomogeneous system, such as partially ordered alloys 
or surfaces with an inhomogeneous concentration profile. 

A more general solution can be found, however, in which the electrostatic potential is 
modified without making effective media for each alloy component in contradiction to the 
assumptions of the CPA. The way to do this is to introduce an additional shift of the one- 
electron potential due to the electrostatic interaction of the electrons inside each atomic 
sphere with the missing charge distributed outside of the sphere and postulate that the 
interaction comes from the boundary between the atomic sphere and the effective medium. 
Such a shift may be associated with an intrasite interaction, which has no connection, at 
all, to the effective medium. 



This is exactly what is done in the locally self-consistent Green's function method [14 



where one goes beyond the single-site approximation for Poisson's equation by means of a 
supercell which models the spatial distribution of the atoms in a random alloy while a CPA 
effective medium is used in the electronic structure calculations beyond a local interaction 
zone (LIZ). If the LIZ consists of only one atom, the LSGF method becomes equivalent 
to the CPA method with a properly defined electrostatic potential and energy In 
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this case, however, each atom in the supercell has its own electrostatic shift given by the 
Madelung potential from all the other atoms in the supercell while the effective medium 
is the same for all atoms. It is clear that such an additional shift for each alloy component 
does not interfere with the CPA because the CPA effective medium is determined on the 
basis of the one-electron potentials including these shifts and because the CPA itself does 
not impose any restriction on the one-electron potentials of the alloy components. 

Following the above arguments two groups have proposed an ad hoc expression for 
the electrostatic shift of the one-electron potential due to non-zero net charges in the 
atomic spheres of the alloy |II^|TT||[2|| . Although the basic models are seemingly different 
and based on different observations, either i), the net charge of an impurity in a metal 
is screened beyond the first coordination shell [16|,[17]] or ii), the net charge of an alloy 



component is proportional to the number of the nearest neighbors of the opposite type 
||10|| , they lead to exactly the same expression for the one-electron potential, i.e., 



v, = (1) 



where Vi is the additional electrostatic shift of the one-electron potential of the 2-th alloy 
component of net charge qi and R\ the radius of the first coordination shell. 

In fact, the models described above are practically identical to the model proposed 
more than three decades ago [I8||l9|] to account for charge transfer effects in the self- 



consistent Hartree scheme based on the tight-binding CPA. In this scheme the variation 
of the z-th atom energy level, <5ej, is proportional to the corresponding charge transfer g,, 
i.e., 5ei = Iqi, where / is some average intra-atomic Coulomb interaction. The non-self- 
consistent limit corresponds to / = 0, while I = oo provides local neutrality ||19|| . In the 
present context one may identify 5ei with Vi and it therefore follows that —e 2 /Ri may be 
considered an intrasite Coulomb interaction. 

Although there is at least some consensus concerning the definition of the additional 
electrostatic shift ([I]), which gives charge transfers quite close to the values obtained in 
supercell calculations [ 20f| , different workers do not agree on the corresponding electro- 



static contribution to the total energy of the random alloy. Some completely deny even 
the possibility of having such a term in a "consistent" SS-DFT-CPA theory pl|,22| while 



others argue about the details of how this term should be defined pO|JT^j23|j2^] . It would 
seem that the presently suggested models of charge transfer effects in the single-site ap- 
proximation to the electrostatic problem, except the trivial elimination of the net charges 
by adjusting the radii of the atomic spheres of the alloy components, may be considered 
neither exact nor even " a consistent theory" . It is the main purpose of the present paper 
to shown that a consistent SS-DFT-CPA theory including a correct description of the 
charge transfer effects does indeed exist. 

Here, we define the electrostatic shift of the one-electron potential and the correspond- 
ing contribution to the total energy in a form which is very similar to that proposed by 
Korzhavyi et al. JT5| , [TTH as well as Johnson and Pinski [|T^], and which provides a prac- 



tically exact solution to the electrostatic problem in the single-site approximation. The 
actual expression for the electrostatic shift in the single-site model for Poisson's equation 
includes one adjustable parameter the value of which is to be obtained in supercell cal- 
culations by the locally self-consistent Green's function method in which the Madelung 
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problem is solved exactly. It turns out that the value of the adjustable parameter is prac- 
tically independent of lattice structure, volume, and alloy composition due to the fact 
that the screening of the electrostatic part of the problem, in the ASA at least, is almost 
universal. 

The paper is organized as follows. In Sec. II we outline the main concepts behind 
the LSGF used in this work and the details of the calculations. In Sec. Ill we present a 
pragmatic solution to the problem of finding the Madelung shift in SS-DFT calculations 
on the basis of the average values of the net charges and Coulomb shifts, < qi > and 
< Vi >, from supercell calculations. We also demonstrate that the linear relation between 
the net charges, q iy and the corresponding Coulomb shifts, Vi, of the alloy components 
discovered for metallic alloys by Faulkner et al. [25] is practically universal in the effective 
medium approach for the Green's function. This means that the response of the electron 
system to the Coulomb field is linear and universal in such systems, and that the screening 
must be universal too. That this is indeed the case is demonstrated in the next section 
where we calculate the distribution of the screening charge in several systems and show 
that it is almost independent of the crystal structure, the alloy constituents, and the 
composition. 

In Sec. V we present a formalism for the screened Coulomb interactions in the single- 
site mean-field approach for the electrostatic potential and energy and demonstrate that 
the conventional assumption of a vanishing Madelung potential and energy is not valid 
in general. Instead, one must include an additional term due to intrasite interactions 
which are, in fact, exactly the screened Coulomb interactions. We also discuss the or- 
dering contribution to the Madelung energy and show why the screening contribution 
may be obtained in supercell calculations for ordered structures. The contribution from 
the screened Coulomb interaction to the generalized perturbation method and the 
formalism is also determined. Finally, in Sec. VI, we demonstrate that the total energy 
of a random alloy may be reproduced exactly in single-site CPA-DFT calculations with 
corrections due to the screening intrasite interaction. 



II. METHODOLOGY 

A. Spatial ergodicity and cluster expansion 

In this paper we will consider only such alloy systems that on an underlying crystal 
lattice with perfect translational symmetry satisfy two conditions: i) spatial homogeneity 
and ii) no correlations between the one-electron potentials at sufficiently large distances. 
For the Coulomb interactions in a random alloy both conditions may be formulated ex- 
plicitly in terms of the average monopole electrostatic potential Vi in the atomic sphere 
around site i due to the charge distribution in the whole of the remaining system. In 
an ordered alloy this potential is the Madelung potential. Specifically, the first condition 
means that any average values of products of potentials must be translationally invariant, 
i.e., 

< VV 3 ...V k >=< T a (V^ ...V k )>, (2) 
where T a is the translation operator T a f(r) = f(r + a), and the second condition is 
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< ViVj . . . VkTaiViVm ...V n ) >=< ViVj ...V k xV l V m ...V n >, (3) 

for a — ► oo. 

According to Lifshitz et al. [2fJ the space formed by the complete set of distinct 



realizations of the potential Vi on the lattice, the operator T a , and the property of spatial 
homogeneity plays the same role in the theory of disorder as the phase space, the operator 
of dynamical evolution, and Liouville's theorem do in statistical mechanics. Moreover, 
according to Birkhoff's ergodic theorem, for any functional /[Vi], where Vi is some random 
realization of the potential on the lattice we have 

lim I / /pT.V;] =< f[Vi\ >, (4) 

i.e., the spatial and the phase-space average are equivalent [p6[l . This equation constitutes 
the principle of spatial ergodicity, according to which all possible finite atomic arrange- 
ments may be realized in a single infinite sample if the conditions (0) and (|j) are satisfied. 

What makes the above principle work in practice is the fact that for self-averaging or 
"measurable" quantities, which have a well-defined limit when the volume of the systems 
approaches infinity, all the correlations of the atomic distribution become unimportant 
at some distances and, hence, the sample may be chosen finite. This may be formulated 



explicitly by means of the cluster expansion theorem |27j which defines the corresponding 



measurable quantity in terms of the site occupation correlation functions 

n = 7r„ + 5>/&> (5) 
/ 

where n q are the coefficients or interaction parameters, £/ =< 5ci5j . . . 5ck > the correla- 
tion function of the figure or cluster / which corresponds to a specific position of the sites 
i, j, and k in the lattice, and Sci = Ci— < q > is the fluctuation of the site occupation 
numbers q taking on values and 1 depending on whether site % is occupied by one or 
the other component. 

According to (|^) there are two practical ways of calculating the properties of a random 
system for which we have = and, thus, H ra nd = i) the cluster or supercell 
approach, where £/ = is satisfied on average only for those clusters / for which 7rj ^ 0, 
or ii) the effective medium approach, which directly gives ttq from some knowledge of the 
alloy components. The first approach is realized, for instance, in the so-called special 
quasirandom structure (SQS) method |28| while the second approach is realized by the 



coherent potential approximation where the real atoms are substituted by a specifically 
chosen effective medium on the lattice. 



B. The LSGF method: a combined supercell-effective medium approach 

The supercell and the effective medium approaches are combined into a single com- 



putational scheme in the locally self-consistent Green's function (LSGF) method [ 14 . 29f] . 
In the LSGF method the supercell approach is used to provide the correct solution to the 
Madelung problem for a given alloy modeled by an appropriate supercell. It is also used 
in part in the electronic structure calculations, which are performed separately for every 
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atom in the supercell by means of the local interaction zone centered at each atom. Inside 
the LIZ the multipole scattering equations are solved exactly, while the region outside the 
LIZ is represented by the effective medium, which is usually taken to be the CPA effective 
medium built on all the one-electron potentials in the supercell. This means that every 
atom of the supercell "sees" only the CPA effective medium outside the LIZ, which ac- 
cording to the CPA definition represents a random alloy. In other words, the one-electron 
Green's function of the supercell (sc) obtained in LSGF calculations may be presented as 



14,29 



G SC = £[G° + £AG^ L / Z L (6) 

i f 

where is the Green' function of the ith atom in the supercell embedded in the CPA 
effective medium, AG if the contribution to G° due to the presence of the specific atomic 
arrangements on the figure / in the LIZ as specified by the correlation function £,^f IZ , 
which is equal to that of the supercell £?? if the figure / is circumscribed by the LIZ, 
i.e., it can be put inside the LIZ in such a way that one of its vertices coincides with the 
central atom of the LIZ, otherwise £k ■— 0. 

Equation (||) clearly shows how the LSGF works, and in particular, how the effective 
medium approach, represented by G®, is combined with the cluster or supercell approach, 
represented by by the second term. It follows from (|6[) that the LIZ allows one effectively 
to cut off the contributions from the clusters which are not circumscribed by the LIZ 
(an equivalent formulation in terms of effective interactions is given in Ref. f29|). If the 
LIZ is single-site, i.e., it consists of only one atom (LIZ=1), the contributions from the 
second term in (Q) vanish and we are left with the usual CPA or pure effective medium 
approach to the electronic structure problem. It is this single-site approach, referred to 
as SS-LSGF, which will be used in most of the present work. It has the advantage over 
the usual SS-DFT-CPA method that Poisson's equation is solved exactly within a given 
approximation for the form of the electron density. 

Note, however, that the LSGF method by no means is restricted to the single-site 
approximation. In fact, it allows us to include local environment effects in the electronic 
structure calculations for the figures circumscribed by the LIZ if on average £,f LIZ> = 
£f C = which is the case in a random alloy. In this respect the LSGF method may be 
considered a self- consistent embedded cluster method (ECM) of the kind proposed by 



Gonis et al. [ 30 1 more than two decades ago. With a proper choice of the supercell used 



to model a given random alloy [59] the LSGF solves two major problems of the ECM: i) 
it provides a set of clusters to represent an alloy with a given short range order and ii) it 
allows one to close the DFT loop with the correct treatment of the electrostatics. 

There is one important point concerning the electronic structure obtained in the LSGF 
with the CPA effective medium which should be mentioned: Despite the fact that it can 
be quantitatively accurate, it is qualitatively different from the electronic structure which 
would result from direct supercell calculations with periodic boundary conditions. That 
is, the electronic spectrum in the LSGF-CPA method is always complex, unless all the 
atoms are equivalent in the supercell (pure metal) or the size of the LIZ is infinite. Thus, 
in the LSGF-CPA the electronic structure of an ordered alloy is never correct, although 
it may be calculated with arbitrarily high accuracy. On the other hand, since Blochs' 
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theorem is not applied to the supercell during the electronic structure calculations, the 
LSGF method is a perfect tool for calculating SQS as opposed to the ordinary band 
structure methods, which in this case lead to a real, i.e., qualitatively incorrect, electronic 
structure of the random alloys due to the fact they present purely a supercell approach. 



C. The choice of the supercell in the LSGF and details of the calculations 

It is possible to obtain the Madelung potential and energy by a combined supercell 
(cluster) and effective medium approach similar to that used in the Green's function 
approach to the electronic structure problem within the LSGF method. However, this 
requires some knowledge of the charge-charge correlations or the screening in the alloy. 
Hence, the only way to solve the problem is to use a supercell model with the Madelung 
potential and energy determined exactly from the bare electrostatic interactions, as it is 
usually done. Here, another problem arises: The supercell should be constructed such as 
to provide zero correlation functions up to the distance where the net charges of the alloy 
components become uncorrelated or completely screened. 

In the calculations presented below we assume at the outset the existence of a short- 
range screening which occurs over the distance of the first several coordination shells. 
This assumption is based on results obtained by the charge-correlated model |ID| , on 
single-impurity calculations |16|JT7|| as well as on the most recent LSGF calculations by 



Ujfalussy et al. |22|| . The latter authors demonstrated that a 16-atom supercell for an fee 
equiatomic random alloy, in which the SRO parameter at the 8'th coordination shell must 
be equal to 1 due to the translation symmetry, i.e., all the atoms in the 8'th coordination 
shell are the same as that at the central site, yields practically the same average charge 
transfer and total energy as a 250-atom supercell, in which the SRO parameter at the 
8'th coordination shell corresponding to a random number generator distribution of the 
alloy components on the lattice should be << 1 (see also fl3l| ). 



Thus, in all the random alloys considered below the distribution of the atoms in the 
supercell was chosen such that the SRO parameters (or pair correlation functions, t^p) 
were exactly zero at least in the first 6 coordination shells and small (not greater 0.01 
in absolute value) up to the 8th coordination shell. Although the multisite correlation 
functions have not been optimized, they should not play a significant role in the ASA 
where only monopole intersite Coulomb interactions are taken into account. 

The electronic structure calculations were in all cases performed in the scalar rela- 
tivistic approximation by the KKR-ASA technique p2| with an s-, p-, and d— basis in 



the framework of either the usual SS-DFT-CPA or the LSGF methods with a CPA effec- 
tive medium. The ASA (no multipole corrections to the Madelung potential and energy) 
has been used in the electrostatic part of the problem. The integration of the Green's 
function over energy was performed in the complex plane over 16-20 energy points on a 
semicircular contour. The local density approximation was used in the DFT part with 
the Perdew and Zunger |33| parameterization of the results by Ceperly and Alder ||34j| . 
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III. NET CHARGE AND MADELUNG POTENTIAL IN METALLIC ALLOYS 



Here we discuss a pragmatic solution to the following problem: Can one devise a 
Madelung potential for the alloy components to be used in SS-DFT-CPA calculations 
such that the charge transfer effects, i.e., the net charges of the alloy components, are 
consistent with those obtained in SS-LSGF calculations where charge transfer effects are 
treated properly? The fact, that such a potential can be found, may seem surprising 
in view of the principle difference between the LSGF and the SS-DFT-CPA methods. 
In the LSGF approach all the atoms in the supercell are different due to their different 
local environment while in the usual SS-DFT-CPA approach one deals only with average 



quantities, i.e., in terminology of the Ref. |3q| , the LSGF supercell approach is equivalent 



to the polymorphous model of the alloy while the effective-medium approach is equivalent 
to the isomorphous model. However, it is obvious, that this can be done on average. 

It was discovered by Faulkner et al. from supercell calculations that the net charges 
on different sites i, qi, and the corresponding Madelung potentials, Vi, obey a linear 
relationship. In Fig. [I] we show such a qV relation for a 512-atom supercell which models 
a random AlsoLiso alloy on an underlying fee lattice. For comparison we also show the 
charge and Madelung shift for AlLi in the ordered Llo structure. All results are obtained 
by LSGF calculations with the CPA effective medium in i) the single-site approximation 
for the electronic part of the problem, SS-LSGF, i.e., LIZ=1, (upper panel) and ii) with the 
perturbation in the electronic structure caused by the local environment up to the second 
coordination shell, i.e., embedded cluster(EC)-LSGF, LIZ=3, (lower panel), included in 
the Green's function. For random alloys the inclusion of more distant coordination shells 
do not affect the results significantly and thus the LIZ=3 results may be considered to be 
converged in the LIZ size. 

The most striking feature of the qV relation obtained in the SS-LSGF calculations 
is the perfect alignment of the qV points along two almost straight lines, one for each 



alloy component. This is, in fact, very similar to what has been observed by Pinski [^G 
in model calculations using the Thomas-Fermi approximation. Furthermore, a change 
of the ratio of the atomic sphere radii of the alloy components, r = Sai/S^, leads to a 
rescaling of the qV points. Hence, for a specific ratio, r = 1.12 in the present case, the 
qV relation collapses into the single point: (q, V) = (0, 0). The existence of this point in 
the SS-LSGF is a consequence of the fact that all Al atoms as well as all Li atoms become 
indistinguishable if the net charges of the alloy components are zero: The difference 
between the atoms caused by local environment effects is solely due to the Madelung 
shift, which is zero in this case. Thus, for this particular choice of r the polymorphous 
model is identical to the usual isomorphous model, or the SS-LSGF method is identical 
to the SS-DFT-CPA method. 

On the other hand, it is clear that the two models are not equivalent when local 
environment effects are included in the electronic structure part of the LSGF method, i.e, 
for LIZ>1. This is demonstrated in the lower panel of Fig. [T] where the local environment 
effects are clearly seen to destroy the strict alignment of the qV points and, as consequence, 
the possibility of choosing electroneutral atomic spheres by a single r value. However, even 
if this were possible, all the atoms, or the corresponding one-electron potentials, would 
still be different. 
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The discussion of local environment effects is beyond the scope of the present paper, 
and the results are included only to demonstrate the qualitative difference between the 
correct results and those obtained by the LSGF method in the single-site approximation: 
Inclusion of intersite correlations in the electronic structure calculations leads to a real 
polymorphous description of random alloys which cannot be mimicked by a single-site 
LIZ. As a consequence, as we will show later the SS-DFT-CPA method can reproduce the 
results of the SS-LSGF exactly, but will, in general, reproduce only approximately the 
correct solution to the supercell or polymorphous model of a random alloy. 

The two qV points for the ordered Ll structure are seen to fall on the qV lines for 
the random alloy as already noted in Ref. P5fl), and, in fact, all the points on the qV 
relation obtained by the SS-LSGF method may be reproduced by a series of ordinary 
SS-DFT-CPA calculations, by using the shift of the one-electron potential defined in a 
way similar to ([[J), i.e., 

Vi = -a^, (7) 

where qi is the net charge of the alloy components, S the Wigner-Seitz radius, and a a 
parameter which may be varied arbitrarily in the SS-DFT-CPA calculations without spec- 
ifying its physical meaning. However, it is important to note, that a = — oo corresponds 
to the electroneutral case (g« = 0) and a = to the limit where there is no response of 
the system to charge transfer effects. As we will see later, the values of the net charges 
qoi obtained in the SS-DFT-CPA calculations with a = are important scaling parame- 



ters. It is also useful to note, that for the Ll structure «li o =0. 8811575 fl37j , and in the 
screened impurity model ([I]) asiM = 0.552669 and 0.568542 for the fee and bec crystal 
structures, respectively. 

Fig. [I] shows the qV relation, indicated by the black line, obtained in the SS-DFT- 
CPA calculations including ([?]) with a varying from —1.5 to 5 together with the SS-LSGF 
results, gray circles. It is clearly seen that the "isomorphous" and polymorphous qV 
relations coincide, and this allows one to make an isomorphous model consistent with the 
polymorphous results. The point is that all the net charges and corresponding Coulomb 
shifts in the polymorphous model have significance only in terms of the average values they 
produce. This is so, because every supercell has its own set of net charges and Madelung 
shifts and, in the case of an infinite system, there is an infinite number of different qV 
points. Their average values, < qi > and < Vi >, however, have a well-defined physical 
meaning as conditional averages of self-averaging quantities, and thus it is the average 
< qV > point which must be reproduced by the isomorphous model. Hence, for a random 
alloy a is given by 

S <Vj> 

errand — ~ • {?) 

e 2 <qi> 

Note, that in a binary AB alloy, it clearly does not matter, for which alloy component, 
i = A, B, a ran d is determined, since < Va > / < qA >=< Vb > / < 1b > P^| - The same 
is true for multicomponent alloys, but in this case, rather than being a trivial consequence 
of the charge neutrality condition, it follows from the physical origin of a ran d, which will 
be discussed in the next two sections. For an fee Al 50 Li5o alloy at S=2.954 a.u. we find 
from the LSGF calculations a ran d =0.60716. 
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Of course, the coincidence of (< qi >,< Vi >) is a necessary, but not sufficient 
condition for the equivalence of the isomorphous and polymorphous models. The two 
models may be called equivalent only if the electronic structure of the random alloy and 
its conditional averages agree. In Fig. || we show that this is indeed the case: The local 
densities of states (DOS) for the Al and Li atoms in AlsoLiso calculated by the SS-DFT- 
CPA method coincide with the corresponding conditional average state densities obtained 
in the SS-LSGF calculations for the 512-atom supercell. For comparison we also show the 
DOS obtained with a =0 corresponding to the "conventional" CPA. Although the latter 
differs from the correct state density, it is obvious, that the neglect of the electrostatic 
shift (|7|) has only a minor effect on the DOS. 

The reason why the average state densities coincide is the following: In the SS-LSGF 
method the difference between the atoms of the same type comes only through the cor- 
responding Madelung shift. A shift in potential leads to a change in the charge transfer 
through a skewing of the local DOS as seen in fig. Therefore, when the conditionally 
averaged DOS is obtained, the skewing contributions from the individual atomic sites 
caused by Vi will cancel and leave only the DOS given by the average < Vi >. Of course, 
this is true only in the SS-LSGF method (LIZ=1). In fact, the local environment effects 
in concentrated random alloys may influence quite strongly the electronic structure of the 
central site of the LIZ. 

To investigate how the qV relation depends on the system we show in Fig. ^ qV 
relations for five different systems including a Cu impurity in Pt (S = 3 a.u.) and four 
random alloys: fee CusoPtso (S = 3 a.u.), fee Al 50 Li 5 o (S = 2.954 a.u.), bec Cu5oZn 50 
(S = 2.7 a.u.), and ternary fee Cu5oNi 25 Zn 5 o (S = 2.65 a.u.). In the plot all charges have 
been normalized by q^ obtained in the no response limit, i.e., a = or Vi = 0, and all 
Madelung shifts have been normalized by qoi/S. To partly simplify the plot we have used 
|goi| m the normalizations, thereby separating the qV relations into two lines rather than 
one. 

The results presented in Fig. ^ show the existence of a universal qV relation. Or 
in other words, the linear-response function \ which gives the change in the net charge 
relative to q 0i caused by Vi, i.e., g« — q 0i = l/e 2 xViS, is a universal constant in metallic 
alloys in the ASA. From the results presented in the figure we find that \ « —0.63. This 
unavoidably leads to the existence of a single, unique a ran d as witnessed by the coincidence 
of all the < qV > points in Fig. [|| Strictly speaking, the slopes of the qV lines are not 
exactly identical and, in fact, a ranc t varies from 0.6 in Cu 5 oZn 5 o and Cu5oNi 2 5Zn 50 to 
0.615 in LiMg alloys, not included in the figure. However, for most practical purposes the 
choice a ran d = 0.607 provides a sufficiently accurate description of the electronic structure 
of random alloys in the SS-DFT-CPA method in the ASA for the electrostatic part. 

IV. SCREENING CHARGE IN METALLIC ALLOYS 

In the previous section we have, in effect, defined a procedure whereby SS-DFT-CPA 
calculations may provide the exact solution to the electrostatic problem in random alloys. 
The only requirements are that the Madelung shift (|7]) is included and that the constant 
errand is obtained from (||) with the average Madelung potential and net charges of the 
alloy components determined in supercell calculations by the SS-LSGF method. In the 
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derivation of the procedure we have used some general arguments which do not clarify 
the physical origin of the universal value of a ran d- However, it is clear, that (0) accounts 
for the missing charge — in the single-site Poisson equation for the i'th atomic sphere. 
Thus, a ran d must be connected to the screening. 

The linear character of the qV relation indicates that the screening in the impurity 
case as well as in the case of a random alloy may be very well described by linear response 
theory. Owing to enhanced electron scattering at opposite regions of the Fermi surface 
linear response predicts in the case of a free-electron gas the existence of long-range Friedel 
oscillations, which however decrease relatively fast (~ r~ 3 ) with the distance. In a random 
alloy, on the other hand, the screening is much more efficient due to the finite life time 
of the Bloch states for the underlying crystal lattice and the spatial distribution of the 
screening density decays exponentially. In this respect charge correlated model [ |10| . |3T|j 
adopted by Johnson and Pinski |12| in the cc- and scr-CPA method or the equivalent 
screened impurity model OJ23Lp|| may be viewed as the first approximations for the 



screening. 

Based on the fact, that a single impurity in a metallic host is a particular case of 
a dilute random alloy, one would expect, and the results for a single Cu impurity in Pt 
presented in the previous section unambiguously indicate this, that the screening effects in 



the two cases are similar. It is therefore surprising that Faulkner et al. [^j and Ujfalussy 
et al. [22| claim that the screening in a random alloy is qualitatively different from that 
found in a single impurity system. In fact, these authors found extremely long-range 
correlations between the Madelung potential at some particular site and the net charges 
at the other sites. Since the Madelung shift on a site is proportional to the net charge on 
the site, this may happen only if there are extremely long-ranged correlations between net 
charges or, in other words, there is no screening. However, this result has been obtained 
on the basis of summations of the direct or bare Coulomb interactions which is, at best, 
an ill-defined procedure, even mathematically. 

To clarify the issue of screening we will perform the following computer experiment 
which will allow us to establish the range of the net-charge correlations or the screening 
in random alloys for one particular site. We set up a 512-atom supercell which represents 
an fee Cu 5 oPt 5 o random alloy (all SRO parameters are equal to zero up to the sixth 
coordination shell and ~ for at least the next 10 coordination shells) and perform self- 
consistent SS-LSGF (LIZ=1) and EC-LSGF (LIZ=3) calculations. We then substitute 
one Pt atom with one Cu atom in some site which, in general, may be chosen arbitrarily. 
However, to keep the atomic distribution as close as possible to the random distribution 
we chose a site the local environment of which corresponds to the random alloy (having 
equal number of Cu and Pt atoms) for the first three coordination shells. We then repeat 
the self-consistent LSGF calculations for the supercell with the substituted atom and 
find new values for the net charges in the supercell. It is clear that, when LIZ=1, the 
difference between the net charges in the two calculations, Ag^, gives the charges induced 
by the change of the net charge at the substitution site (in the case of LIZ=3 the local 
environment effects also effect the charge transfer). This charge is simply the screening 
charge. 

In the upper panel of Fig. |5| we have plotted the normalized, induced charges 
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at the first eight coordination shells around the substitution site % — for Cu-Pt substi- 
tution in a CusoPtso random alloy and for Cu-Pt substitution in pure Pt (S — 3 a.u.) 
obtained in the single-site approximation for the electronic structure (LIZ=1) as well 
as with local environment effects included (LIZ=3). One may see that, while the local 
environment in a pure metal hardly affects the distribution of the net charges, it does 
introduce a dispersion in the distribution of the net charges in the random alloy, which is 
quite substantial at the first coordination shell but which practically disappears beyond 
the 5'th coordination shell where, in fact, all the induced charges almost vanish. 

To demonstrate that the net charge of the Cu " impurity" indeed becomes screened we 
show in the lower panel of Fig. [| the total normalized induced charge in the z'th shell 

QL = X>Qi, (io) 

j=0 

where Zj is the coordination number of the j'th shell. I is seen that Qi vanishes beyond 
the 7'th coordination shell in all cases, and we conclude that the screening in a random 
alloy in the single-site approximation is practically the same as the screening in the case of 
a single impurity in a pure metal. There is neither qualitative nor quantitative differences 
between the impurity and the alloy cases. 

In the upper panel of Fig. ^ we show the distribution of the screening charge (not to be 
confused with the screening density: the screening charge is, in fact, the screening density 
integrated in the corresponding atomic sphere) for a Cu impurity in fee, bec, and bet Pt 
plotted as a function of the distance from the impurity site in units of the Wigner-Seits 
radius, S. It is clear that the screening charge follows a single, common curve which does 
not depend on the structure. In fact, by changing the c/a ratio in the bet structure one 
may completely fill the remaining gaps in the calculated curve. In the lower panel of the 
figure we have collected the results for the distribution of the screening charges in seven 
different systems including such hosts as Pt, Al, Cu, V, Na, and K. It appears that the 
screening in metallic alloys depends neither on the crystal structure nor on the nature of 
the alloy components, at least, when described within the ASA. 

The universal picture of the screening in alloy systems found above is partly destroyed 
when the electrostatics is treated more correctly, for instance by including multipole 
moment contributions to the one electron Madelung potential and energy. However, the 
ASA still gives a qualitatively correct picture and catches the main physics behind the 
phenomenon. Hence, it is worth to discuss the origin of such a universality in both the 
screening and the response function. 

First of all, it was understood long time ago that the net charges in the atomic spheres 
of the alloy components has very little in common with the "charger transfer" in terms 
of the redistribution of the electron charge between the alloy components (see, for in- 
stance Ref. |38[]). Even in the case of the so-called ionic solids the self-consistent charge 
distribution is very close to that obtained from a linear superposition of the free-atom 
electron-densities |3TI] and this is the reason for the success of the charge-correlated model 



|10| , in which the net charge is proportional to the number of nearest neighbors of the 



opposite type (see also Ref. 
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What we are seeing is basically a size effect: The net charges originate from the 
redistribution of the electron density in the interstitial region between the atomic spheres. 
The electron density in metals and their alloys in this region is very smooth and may be 
well described by a free-electron model, even for transition metal alloys. The interstitial 
density is much easier to perturb than the density closer to the atomic nuclei and it 
participates in the screening. On this basis one may, in fact, developed a model based 
on linear response theory which leads to a semi-analytical description of the universal 
screening. However, this is beyond the scope of the present paper. 

The calculated distribution of the screening net charge may be used to obtain the 
screening contribution to the one-electron potential in the single-site model for the Poisson 
equation given by (0) with a equal to 

e i Ki 

where Ri is the radius of the z'th shell with coordination number Zj. Using the results for 
the screening charge, Q iy in the case of Cu-Pt substitution for one particular site (not on 
average!) in the CusoPtso random alloy we find from (|TT| ) after summation up to the 8'th 
coordination shell that a scr = 0.60572. At the same time, the average values of the net 
charges and the Madelung potentials, < q > and < V >, in conjunction with (H) gives 
errand = 0.60530. That is, a ran( i = a scr . Since, there is only one effective medium in the 
supercell LSGF calculations, it is obvious that the screening is the same for all the sites 
and, thus, should not depend on the alloy component, that is, the ratio < Vj > / < qi > 
does not depend on the alloy components. 

It is important to understand that a ran d appears in the formalism due to the intrasite 
interactions between the electron density inside an atomic sphere and its screening (or 
missing) charge. Therefore the inclusion of (]?]) in SS-DFT-CPA calculations does not in 
any way contradict the "conventional" mean-field picture according to which the contri- 
bution from intersite interactions to the Madelung potential is zero in the SS-DFT-CPA. 
In fact, neglecting the term will lead to incorrect results. On the other hand, the param- 
eter cxli which allows one to perform SS-DFT-CPA calculations identical to those of the 
SS-DFT-CPA for the Ll structure, has an origin different from that of a scr determined 
by flllD . The former is a constant, which apart from the Ll symmetry does not depend 
on anything, and in particular, not on the screening, while the latter is a system depen- 
dent parameter which is entirely determined by the type of intersite interactions or the 
screening in a given system. 

This means that there is no connection between the Madelung constants for ordered 
structures and a scr . The reason why it is possible to obtain ot scr from the supercell 
calculations for the completely ordered structures, as we have just done, is the fact that 
at the large distances, where the atomic-distribution correlation-functions are not zero 
anymore, the "real" intersite interactions do not contribute to the electrostatics due to 
the short-ranged screening. We will return to this point in the discussion of the Coulomb 
energy of a random alloy, but here we would like to comment on the use of the single-site 
approximation in the Green's function calculations. 

Our LSGF calculations of impurities in different metals indicate that the problems 
observed in the single-site Green's function impurity calculations by Stefanou [17]] and 
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Drittler et al. j4(| do not originate from the single-site approximation for the Dyson 
equation, but from Poisson's equation, which these authors also solve in the single-site 
approximation. The use of the screening electrostatic shift for the one-electron potential 
allows one to solve the impurity problem in the single-site Green's function formalism in 
the ASA or in the MT approximation almost exactly. This is so because the impurity case 
corresponds to the dilute limit of a random alloy where the concentration of one of the 
alloy components approaches zero. In this case, the contribution to the electronic structure 
due to the local environment effects becomes negligible, and the electronic structure of 
the impurity obtained by the single-site or the cluster Dyson equations become almost 
identical. Such an effect may be seen, for instance, in Fig. ^ The dispersion of the 
screening net charge found in the concentrated alloy case, i.e., CusoPtso (LIZ=3), vanishes 
in the case of a Cu impurity in Pt. The effect is, in fact, the origin of the increasing 
accuracy of the CPA with decreasing concentration of one of the alloy components. 

V. SCREENED COULOMB INTERACTIONS 

A. The Madelung energy of a random alloy in the single-site mean-field 

approximation 

The formalism describing the electrostatics of random alloys in the single-site mean- 
field approximation, where all the A and all the B atoms are represented only by the 
appropriate conditional averages, is indeed a trivial one, although it has remained quite 
a confusing issue for more than 30 years, starting from the definition of the electrostatic 
energy of a random alloy given by Harrison in 1966 in connection with pseudopotential 
theory ||41|| . According to Harrison the Madelung energy, which may be associated with 



the electrostatic interactions of the net charges of the ions of a random A C B!_ C alloy, is 
given by |^1 



E%t" = -yf c(l " c)(Z A - Z B f + e ^f~Z\ (12) 

where Zi are the ion charges of the alloy components, ocm is the Madelung constant of the 
underlying lattice, S the radius of the Wigner-Seitz sphere, Z the average charge equal 
to cZa + (1 — c)Z B , and a some constant. 

Equation (|TjJ) was, in fact, not derived specifically for Coulomb interactions. Rather, 
it was suggested on the basis of a more general consideration of the band-structure con- 
tribution to the total energy of a random alloy within second order perturbation theory. 



Later, it was shown by Krasko [52] that in a system with randomly distributed A and B 



ions, having charges Z4 and Z B , respectively, embedded in a medium of compensating 



charges, the electrostatic energy included only the second term in ( |T2| ) while the first term 
vanished. 



It is obvious that (|1^) is valid also in the CPA-DFT if one substitutes the ion charges Zi 
by the net charges qt of the atomic spheres. In the ASA the average charge cqA + (1 — c)q B 
is zero, and in this particular approximation the contribution to the Madelung energy of 
a random alloy given by the second term in (|12"D will also vanish. Note, however, that the 
term vanishes neither in inhomogeneous systems [050] nor i n the MT-approximation || . 
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This result was criticized by Magri et al. [T(| who a decade ago deduced from cal- 
culations for ordered compounds that the net charge of the atomic sphere of an alloy 
component is proportional to the number of nearest neighbors of the opposite type, and 
on this basis developed the charge correlated model in which the Madelung energy of a 
binary random alloy is given by 

ETat SS = -y^ 1 - <0(«A ~ q B f, (13) 

where a cc is a constant equal to 0.5947049 and 0.60817846 for the fee and bec structures, 
respectively. Exactly the same result has been obtained later by Korzhavyi et al. [ |15l , |Tl"ll 



and by Johnson and Pinski [12|] in their "screened" models for the single-site CPA-DFT. 
The difference between these models lies only in the way the parameter a cc is determined 
(from 0.4397212 for bec in Ref. @ to 0.54282038 for fee in Ref. gt|). A discussion of 



the issues involved may be found in Refs. |T2"| , P3"| , P4T , |3~1"| , P0" 



It is certainly surprising, that ([13]) is exactly the first term in (0) which was shown to 
vanish in the case of two charges, here and randomly distributed on an underlying 
lattice in a compensating homogeneous effective medium. This means that either the 
models are inconsistent with general theory or that there should be some reason, not 
accounted for in the derivation by Krasko, for the presence of the first term in (0). 

Assuming the existence of only on-site and pair-wise interactions the Hamiltonian of 
a binary A c Bi_ c alloy may be written 

H = Ei&R + " c fl )] + [vbxCsPx + - c R )c R/ (14) 

R / R^R' 

+ vj^,C R {\ ~ CBf ) + «f£,(l - C fl )(l - CBf )], 

where are on-site or intrasite interactions, which we will assume depend only on the 
type of atom on site R, v^, are pair potentials acting between X and Y atoms at site R 
and R', respectively, and c R is the site-occupation operator taking on the value 1 if there 
is an A atom on site R and otherwise. Using Scr, defined by cr = c + 5cr, we may 
rewrite the Hamiltonian in the equivalent form 

H =\H Vrr'Sc r 8c r , + W [c\ a r b r , + 2c(l - c)v A R B R , + (1 - c) 2 «], (15) 

Z R,R' Z R+R! 

where the first term includes the intrasite interaction [R = R') 

VR=o = 2[-^—e A + -e B } (16) 
(1 - c) c 

as well as the intersite interactions (R ^ R') 

Vrri=v aa +v b r b ri -2v ab ri . (17) 
Upon Fourier transformation of the first term we find 

H=-g- I dq V(q)c q c* q + ~ £ [c 2 v AB + 2c(l - c)v£ B + (1 - cfv BB ], (18) 
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where the second term is the average contribution to the energy due to pair interactions 
which in the case of direct Coulomb ion- ion interactions is q x q y /\H— R'| combined with the 
corresponding contribution from the interaction between the ions and the homogeneous 
compensating charge. This is exactly the second term in (]T2|) . 

The first term in (|IBD is usually associated with the configurational contribution to the 
energy of the system, but this is correct only if the contribution from intrasite interactions 
is zero. It is easily evaluated in a completely random alloy, where all the occupation num- 
bers are uncorrelated and therefore c q c* = c(l — c)/N (which provides the normalization 
of the ordering energy per atom). One finds 

/ dq V(q)c q c* q = c(l - c) / dq V(q) = n BZ c(l - c)V R=0 , (19) 

JBZ JBZ 



c e 



o 



which according to flTID is equal to ce$ + (1 

It now remains to define the on-site interaction term e l , which results from the in- 
teraction of the net charge g« in the alloy with the corresponding screening charge, in 
such a way that and the corresponding on-site Coulomb potential Vi given by (0) are 
consistent within DFT, i.e., Vi = dey^qf. 



e 2 a 



scr 2 



qf. (20) 



2 S 

Using this definition the first term in the Hamiltonian ( f[8|) may be written 

ETat™ = E Mad SS = -J^flcql + (1 " c)ql] = -^c(l - c)^f(q A - q B f (21) 
= c(l-c)V scr (R = 0), 



which is exactly the first term in ( |T2"D and which, according to the above derivation, arises 
from on-site or intrasite interactions such as the screening interactions in metallic alloys. 

It is now clear why the first term in flT^) is absent in the work by Krasko: His derivation 
is based on inter-site Coulomb interactions only. Thus, Krasko's result is valid in the 
absence of the screening intrasite interactions, which is the case, for instance, in the 
pseudopotential formalism within second order perturbation theory j|T|. However, at the 



same time Harrison's result is in fact correct under the more general assumptions needed 
in the SS-DFT-CPA calculations. 



B. The configurational part of the Madelung energy and potential 

The previous section may seem trivial: First we define our inter- and intrasite interac- 
tions in real space, then we Fourier transform, and from the Fourier transform we return 
to the initially defined intrasite term. However, in some formalisms, such as pseudopo- 
tential theory or the formalism jl^-|55|], V(q) is already defined and this may lead 
to problems with the correct definition of the configurational part of the total energy. 
The point is that the intrasite interactions do not contribute to the configurational part 
of the total energy which in real space may be written as (here, we do not consider the 
contribution from multi-site interactions) (46| 
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Hconf = - v RR'$ c R Sc R'- ( 22 ) 

1 R^R' 

Therefore, if the configurational Hamiltonian is written in terms of V(q), e.g., in the con- 
centration wave formalism, it must be corrected by the subtraction of the corresponding 
intrasite interaction, i.e., 



H, 



If 1 

conf = ^ — / dq V(q)c q c* - -c(l - c)V R=0 (23) 

= — / d q i v (q) - v R=0 ]c q c* 

L\Ibz jbz 

where we have used the sum rule for the concentration wave density c q c*: 
Jbz dq V(q)c q c* q = Vt BZ c(l - c). 

The subtraction of the intrasite term in ( p3"D is crucial for obtaining the correct ordering 
energy in pseudopotential theory and the formalism p3|-^1 as well as for making the 
whole theory consistent. Let us, for instance, consider the Madelung energy of a binary 
completely ordered alloy with two non-equivalent sublattices. It is easy to show that its 
Madelung energy can be presented exactly in the form flL2"|). For instance, the Madelung 
energy of the Ll ordered phase is 

ESS, = -J^ 1 ~ c )(^ - **? + 7^ ( 24 ) 

where the last term is zero in the ASA since q=0 as in the random alloy case, but 
now aLi is a constant which appears due to the intersite Coulomb interactions. In the 
Appendix of Ref. P7J] it is shown that, in fact, aLi (qA — 9b) 2 '/ 'S is the Fourier transform 



of the effective direct electrostatic interaction at the corresponding superstructure vector 
k Llo = 27r/a(100), i.e., V es (k Llo ) = a Llo (q A - q B ) 2 /S. 

On the other hand, ma y a l so be found as the sum of the electrostatic energy 

of the completely random alloy, E r A ^~ ss , given by ( |2"T| ) and the ordering energy, AU: 
E Mad = E Mat ss + Au '■ Since the Madelung energy of the ordered Ll alloy, E^ d is 
uniquely defined in terms of the corresponding Madelung constant, which has nothing to 
do with the screening in the alloy, it is obvious that such a screening term must be present 
in the ordering energy, AU, to compensate the screening contribution. 

Indeed, as shown in the Appendix of Ref. [IB| the ordering energy in the Ll structure 
can be written in the form as 

AU Mad = l v 2 [V es (k Ll0 ) - V scr (R = 0)], (25) 

o 

from which it is easy to see that in the in the completely ordered state, where the long 
range order parameter r]=l, the last term in ( fZSD is exactly the Madelung energy of 
random alloy at the stoichiometric composition (c(l-c) = 1/4), and thus E r ^~ ss + AU = 
l/8a Llo (q A - q B ) 2 /S. 

This illustrates an important point: The ordering energy represented in reciprocal 
space in the concentration wave formalism must be corrected by the exclusion of the 
intrasite term, otherwise the theory will not be consistent. Equation (p3"D gives the correct 
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definition of the ordering energy considered more thoroughly in the Appendix of Ref. |46| . 



The intrasite interaction must also be subtracted when one considers the energy of SRO 
effects, and thus the correct Krivoglaz-Clapp-Moss expression must have V(q) — Vr =0 



instead of V(q), which is exactly the case in Krivoglaz's derivation [j47|l - Note, however 



that this problem does not exist if the Krivoglaz-Clapp-Moss expression is used together 



with the so-called Onsager correction [[Hj] under the condition that it is properly defined. 

The reason, why it was possible to calculate a scr on the basis of the ordered structures, 
is the fact that, in an ordered binary alloy with only two non-equivalent sublattices, 
one has an exact cancellation of the screening contribution to the Madelung energy and 
potential. This does not happen, however, in the general case of a supercell with n > 
2 non-equivalent sublattices. Here, the Madelung energy may written as sum of the 
contribution from the intrasite screening interactions 

E s M~r = E = £ & + (! - c )4- E &] ( 26 ) 

2yv b i 2 b N A i=A Mb i=B 

[^ e ^(c< q A> 2 +(l-c)<q B >% 

where N A and N B are the number of A and B atoms, respectively, and the ordering energy 
due to the intersite interactions 

e 2 

AU Ma d = ^ E 7<(«k - ow)Agf.. (27) 

Here, 7, is a normalizing coefficient, a constant due to the bare electrostatic interactions 
between the net charges for the superstructure vector ki which may be calculated from 
the Madelung constants a M of the corresponding supercell similar to the an considered 
above, and Ag^. the difference between the charges in the crests and in the troughs of the 
concentration wave in the supercell. In the case of a binary alloy with two non-equivalent 
sublattices, there is only one ki and Aq^ = (q A — 

If depends only on the structure and describe the bare electrostatic interaction 
between net charges, then Ag%. "dresses" these interactions according to the real charge 
distribution in the alloy (an equivalent description in real space in the charge correlated 
model is given by Wolverton and Zunger ||8j , who also show that the Madelung energy of 
the random alloy has intrasite character). If the net charges in the supercell are screened 
(or uncorrelated) at distances less than half the period of the concentration wave with 
wave- vector ki then Aq^—0 and the corresponding contribution to the ordering energy 
vanishes. If the supercell includes only long-range concentration waves, the corresponding 
ordering contribution to the Madelung potential and energy becomes zero. 

Let us finally mention the fact that the Madelung energy of a random alloy obtained in 
supercell calculations (p6|) is not equal to the Madelung energy in the single-site calcula- 
tions, and thus it cannot be used to obtain a scr . The reason is simply that the Madelung 
energy is not a self-averaging quantity. However, the Madelung potential is, and it is clear 
that 

< V X >~= E ft* = ^ < ix >= VT, (28) 

6 i=x b 
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which allows one to use (H) to obtained a scr in the supercell calculations and shows why 
errand is exactly equal to a scr . 

C. Intersite screened Coulomb interactions 

Although the screened Coulomb interactions have an intrasite character, they may 
contribute to the effective pair intersite interactions of the kind obtained in the generalized 



perturbation method (GPM) ||T9| , f49l1 because the screening charge is located on several 



of the coordination shells around each atom. This was, in fact, already recognized by 



Ducastelle flnj who derived the contribution to the GPM potentials from the screened 
Coulomb interactions in the framework of the Hartree-Fock tight-binding CPA theory. 

The existence of an additional electrostatic term due to the screening is also consistent 
with Andersen's the force theorem [|50|| , which states that the change in the total energy 
of a system due to some perturbation to first order is given by the change in the sum 
of the one-electron energies obtained from frozen one-electron potentials plus the change 
of the electrostatic energy due to the perturbation. In fact, this latter contribution from 
the screened Coulomb interaction has been completely neglected in a number of first- 
principles calculations of GPM interactions |51|-p4|. Here, we will therefore show how the 
screening contribution to the GPM potentials may be defined and obtained on the basis 
of the calculated spatial distribution of the screening charge. 

GPM-like pair interactions, usually defined by (|PT| ) for a specific lattice vector R, may 
be determined as the site-projected part of the change in the total energy when two atoms 
of different types in a completely random alloy are exchanged between sites infinitely far 
apart in such a way that their neighbors at the relative position R are of the opposite 
type after the exchange. This is schematically illustrated in Fig. [?|. That part of the total 
energy which should be accounted for is half the site-decomposed total energy written in 
terms of the intersite interactions or interatomic potentials, i.e., 

V^\R) = l -[Ef\R)-Ef\R)}. (29) 
Here, E^ is the total energy due to pairwise interactions of the unperturbed system pro- 

(2) 

jected onto site and E 2 is the same quantity after the exchange. A similar expression 
is also valid in the case of multisite interactions, but this will involve a more complex 
exchange of atoms and will not be considered here because the screened Coulomb inter- 
actions do not contribute to the effective multisite interactions in the ASA. 

Within multiple scattering theory as well as in the tight-binding approximation a 
Green's function formulation allows both site- and "path" -decomposition of the electron 
density and thereby makes it possible to write down an analytical expression for the one- 
electron contribution to the n-site interactions, V^ n '(R), in the CPA |l9| , f49"fl . Concerning 
the screened Coulomb interactions one must, however, proceed differently. There are 
several ways to do so, but here we will present a straightforward approach. 

In the sense of the CPA and single-site mean-field theory we will use an effective 
medium approach, assuming that at all sites, i.e., within the atomic spheres assigned to 
each site, there is an electroneutral effective medium except at the two sites and R 
under consideration. In those two sites we must use the actual values of the net charges 
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of the alloy components, which in the effective medium approach, are the average net 
charges, q A and qs, of the alloy components. 

In the first-principles methods, however, these net charges depend on the specific 
choice of the size of the atomic spheres and thus they must, in principle, go together 
with the corresponding screening cloud. Since we calculate the change in the electrostatic 
energy of the two systems shown in Fig. |7] projected onto site due to the exchange of 
A and B atoms in positions R, we must include only the interaction of the net charge at 
site with the net charge at site R and its screening charge. That is, the interaction of 
the net charge with its own screening charge must be excluded as it is included in the 



definition of the screened on-site interactions, see (|TT|). Thus, the first term in ( p9|) for 



the system before the exchange of atoms has been made, E{ cr (R), is 

ET(R) = e 2 q A E ^ + #1* E ^F- (30) 

RVO R'^0 U 

Here, q^i is either the net charge of the ith component, q^ if R' = R or the corresponding 
screening charge if R' ^ R. A similar expression may be written for i?| cr (.R), after the 
exchange of the A and B atoms in the R-sites, i.e., 

ET r (R) = eV E + ^ E ^r- (31) 

RVO n RVO U 

The resulting expression for the screened Coulomb interactions which should be added to 
the usual one-electron term is therefore 

VUR) = f lu £ ~ iBW ~ ib £ qBR, R ," AR ' ] (32) 

/ RVo n'yto n 

e 2 g(|R-R|) 

= 17{ ( Ia-Qb) 2^ ™ ' 

Z R'^0 n 

where Q(R) is the normalized screening charge defined in (||), and where we have used the 
condition that the screening does not depend on the type of the atom. Finally, performing 
the summation in ( p2[ ) one may define the screened Coulomb interactions as 

V SC r(R) = j(qA-qB) 2 ^j^. (33) 



It is easy to see from (|32|) and (|Tl|) that a scr (R = 0) = a scr = a ran d and therefore, 
V scr (R = 0) is exactly the on-site screened interactions that defines the Madelung energy 



of the binary alloy which has exactly the same form (]2T|) . This on-site interaction must 
be included in the definition of the interactions ||43|| , as has been demonstrated in 
the previous section (see also Ref. [0). When R ^ 0, V scr (R) defines the intersite 
screened Coulomb interaction contribution to the GPM-like effective interactions. Since 
the screening in the ASA, is practically universal these interactions have the universal 
form presented in Fig. H. 
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VI. THE TOTAL ENERGY IN THE SINGLE-SITE CPA AND THE 
SUPERCELL LSGF METHODS 



The fact that the Madelung energy of a random alloy described either by the effec- 
tive medium model defined by the SS-DFT-CPA method or by the supercell model in 
conjunction with the SS-LSGF method differ from each other, has neither consequences 
for the final result for the total energy of the random alloy nor even for the partial and 
local contributions to the total energy. This follows simply from the fact that the density 
of states and its average local contributions are the same in the two methods, as shown 
above. 

In Table | we compare the total energy and its components in a CusoPtso random 
alloy calculated by the SS-DFT-CPA method with a scr =0.60572 and by the SS-LSGF 
method on the basis of a 512-atom supercell, in which the atomic positions of Cu and Pt 
have been chosen such that the SRO parameters are equal zero at the first 7 coordinations 
shells (LSGF-1). The agreement between the two calculations is seen to be excellent if 
one combines the electron-nucleus, the electron-electron, and the Madelung contributions 
to form a total Coulomb energy, E coui = E d _ nuc + E e ^ e i + E Mad . 

The accuracy of the SS-DFT-CPA method with the appropriate screening contribu- 
tion to the Madelung potential and energy may be appreciated if one compares the results 
of a 512-atom supercell calculation performed by the SS-LSGF method (LSGF-2) where 
the distribution of the Cu and Pt atoms have not been optimized after the application 
of the random number generator leading to quite small, but not zero, SRO parameters. 
The values of the SRO parameters for the first 7 coordination shells are —0.005208 (1), 
0.026041 (2), 0.007161(3), -0.014323(4), -0.021484(5), 0.0390625(6), -0.0136718 (7), re- 
spectively, which are approximately the same, as in the LSMS calculations in Ref. ||55|| . 
The agreement between SS-DFT-CPA results and SS-LSGF calculations with a properly 
chosen supercell (LSGF-1) is obviously better than between two SS-LSGF calculations. 

VII. CONCLUSION 

The screened Coulomb interactions which are due to the interaction between the net 
charge of an alloy component and its screening charge must be included in a consistent 
single-site mean-field theory of the electrostatics in random alloys. In this paper we 
have shown how this may be done and we have calculated the spatial distribution of the 
screening charge which in the ASA is found to be practically universal for homogeneous 
systems. A formalism that describes the contribution from for the screened Coulomb 
interaction to Madelung potential and energy as well as to the effective interactions of 
the GPM-type is presented. 
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FIG. 1. The distribution of the net charges, qi, and corresponding Madelung potentials, 
Vi in the 512-atom supercell, modeling a random AI50L150 alloy, ordered Llo alloy and in the 
single-site CPA-DFT calculations obtained by varying a. 




FIG. 2. The local density of states in the AhjoLiso obtained by the LSGF method with LIZ=1 
and by the single-site CPA-DFT with different values of otrand- 
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FIG. 5. The distribution of the screening charge in the random CusoPtso alloy and in Pt 
due to substitution of a Cu atom. 
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FIG. 6. The distribution of the screening charge in different metals having different crystal 
structure and lattice parameter. 
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FIG. 7. Two systems, whose O-site projected Coulomb energy is to be used in the calculation 
of the screened effective interactions at distance R. 
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FIG. 8. The intersite screened Coulomb effective interactions obtained from the normalized 
screening charge presented in Fig. 6. 
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TABLES 



Site 



Energy 



ss-CPA-DFT 



LSGF-1 



LSGF-2 



Cu 



Kinetic 

< Eel—nuc ^* 

< E d - e i > 

< Euad > 

< E C oul > 

< E xc > 

< Ecu > 



3360.076110 
-7974.160832 

1439.080777 
-0.004193 
-6535.084248 
-130.026085 
-3305.034222 



3360.076294 
-7974.157257 

1439.078359 
-0.005646 
-6535.084544 
-130.026000 
-3305.034250 



3360.077674 
-7974.178557 

1439.099272 
-0.005994 
-6535.085279 
-130.026621 
-3305.034226 



Pt 



Alloy 



Kinetic 

^ Eel—nuc ^* 

< E e i^ e i > 

< E Mad > 

< Ecoul > 

< E xc > 

< E Pt > 



E, 



tot 



42188.794140 
-92747.101049 

14378.917863 
-0.004193 
-78368.187379 

-693.866856 
-36873.260095 



42188.794273 
-92747.093691 

14378.911671 
-0.005533 
-78368.187553 

-693.866786 
-36873.260068 



-20089.147159 



-20089.147159 



42188.791806 
-92747.030284 

14378.849707 
-0.005472 
-78368.186049 

-693.865919 
-36873.260162 



-20089.147194 



TABLE I. The total energy (in Ry) of CusoPtso random alloy and corresponding con- 
tributions obtained in three different calculations: by the single-site CPA-DFT method 
(ss-CPA-DFT), in the 512-atom supercell LSGF calculations with optimized atomic dis- 
tribution, providing zero SRO parameters up to the 7th coordination shell (LSGF-1), 
and with atomic configuration immediately after random number generator (LSGF-2). 
(Ecoul = E el _ nuc + E e i_ ei + E Mad ) 



28 



